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Abstract 

We investigate the ground state of a uniaxial ferromagnetic plate with perpendicular easy 
axis and subject to an applied magnetic field normal to the plate. Our interest is the asymptotic 
behavior of the energy in macroscopically large samples near the saturation field. We establish 
the scaling of the critical value of the applied field strength below saturation at which the ground 
state changes from the uniform to a branched domain magnetization pattern and the leading 
order scaling behavior of the minimal energy. Furthermore, we derive a reduced sharp-interface 
energy giving the precise asymptotic behavior of the minimal energy in macroscopically large 
plates under a physically reasonable assumption of small deviations of the magnetization from 
the easy axis away from domain walls. On the basis of the reduced energy, and by a formal 
asymptotic analysis near the transition, we derive the precise asymptotic values of the critical 
field strength at which non-trivial minimizers (either local or global) emerge. The non-trivial 
minimal energy scaling is achieved by magnetization patterns consisting of long slender needle- 
like domains of magnetization opposing the applied field. 
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1 Introduction 

Ferromagnetic materials offer a fascinating example of physical systems capable of producing an 



extraordinarily rich variety of spatial patterns 16 . By a pattern in a ferromagnet, one usually 



understands a stable spatial distribution of the magnetization vector in the sample. This definition 
reflects the mesoscopic nature of the magnetization patterns: they are observed on the length scales 
significantly exceeding the atomic scale (making the definition of the magnetization per unit volume 
meaningful), yet they are susceptible to small random fluctuations due to thermal noise, with the 
noise providing a selection mechanism for observable patterns. 

On the mesoscopic level, the theory describing the spatio-temporal dynamics of the magnetization 
patterns in ferromagnetic materials is formulated in terms of partial differential equations (with a 
possible addition of stochastic forcing [2]) for the magnetization vector M = ]V[(r,t) [16,23,25]. At 
the center of the theory is the micromagnetic energy functional £[M.] describing the contributions 
of different physical interactions (for specifics, see the following section) [IT 16 19]. Magnetization 



patterns are viewed as global or, more generally, local minimizers of forming mainly due to 
the competition of the exchange, anisotropy, and the magnetostatic interactions, with the applied 
external field playing a significant role 16p3 . Because of the non-local nature of the magnetostatic 



forces, their effect can depend significantly on the geometry of the ferromagnetic sample 11 12l 16 
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In bulk crystalline materials the local anisotropy energy and the short-ranged exchange energy act 
jointly to favor magnetization distributions in the form of extended magnetic domains in which the 
magnetization vector stays nearly constant, separated by domain walls, where the magnetization 
direction changes abruptly. It was already realized in the pioneering works of Landau and Lifshitz 



24 and Kittcl 18 that, while the structure of the domain walls may not be significantly affected by 
the long-range magnetostatic forces, these forces should determine the relative spatial arrangement 
of the domains with different orientation of the magnetization. In fact, since the total magnetostatic 
energy scales faster than volume as the size of the system increases, in large samples the effect of 
long-range magnetostatic interactions becomes dominant. As a result, the magnetization patterns 
develop rapid oscillations to cancel out the induced magnetic field and form intricate structures, 
which are generally referred to as branched domains, even though the actual topological branching 
of the domains is not really required. 



Despite a long history of observations of branched domain structures in ferromagnetic materials 16 
and related systems (see e.g. 23p9 32]), mathematical understanding of the branching phenomenon 



started to emerge only recently with the ansatz-free analysis of energy minimizing structures [5j[7j 
(there is, of course, an extensive literature of ansatz-based studies, see e.g. 15][l7]l28] ) . In particular 



for bulk crystalline ferromagnets in the absence of an applied field the first rigorous analysis of the 
branched domain structures was performed in the work of Choksi and Kohn [6]. They studied 
a sharp interface version of the micromagnetic energy and were able to obtain matching (in the 
sense of scaling with the sample thickness) upper and lower bounds for the energy of minimizers 
of the reduced energy. We note that the connection of the sharp interface energy to the full 
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micromagnetic energy in the Hmit of high anisotropy was recently estabhshed in 27 . The results 
of are suggestive that the energy minimizers of the sharp interface micromagnetic energy are 
in some sense not very different from the branched domain ansatz used as a trial function in the 
calculation of the upper bound of the energy of the minimizers. The latter shares many common 



features with the branched domain structures observed in experiments 16 . Since then, similar 
results have also been obtained for models describing type-I superconductors in the intermediate 
state [Slls] and diblock-copolymers undergoing microphase separation 26 . 

Note that the presence of a moderate applied magnetic field does not alter the situation qualitatively. 
On the other hand, if a very strong external magnetic field is applied to the sample, then it 
will obviously overwhelm all other effects and result in a uniform magnetization pattern in the 
direction of the applied field. It is then clear that a bifurcation from the uniform to a non-uniform 
magnetization pattern will occur when the field strength is gradually reduced. Let us point out that 
this transition would typically occur via nucleation and growth of new domains and is, therefore, 
accompanied by a hysteresis. In other words, in a certain range of applied fields one should find 
coexistence of different types of patterns. Their relative stability and the transition pathways 
between them are, therefore, important questions to be addressed. Note that these questions also 
naturally arise in various other problems of energy driven pattern formation, such as type-I and 



type-II superconductors and Ginzburg-Landau models with Coulomb repulsion 111 5 8 26 . 



Main results 

We investigate the properties of the magnetization patterns in bulk uniaxial crystalline ferromagnets 
in the presence of external magnetic field applied along the material's easy axis. We are interested in 
the transition to non-trivial energy minimizers occurring near the saturation field in ferromagnetic 
plates with perpendicular easy axis. 

In Section [3j we establish the scaling behavior of the minimal energy in dependence of the plate 
thickness and the exterior field. The precise result is stated in Theorems |3.1| and |3.2[ As a 
consequence, we get that in macroscopically large plates the transition from the monodomain to 
the branched magnetization pattern occurs when the strength Hcxt of the applied field i?oxt satisfies 

Hs-H,.t - TiTTTT^r— — In — — ^ , (1.1) 



L^Ms{K + AnM^) \ AK 
where Hs — AnMs is the saturation field. We refer to Section [2] for the precise definitions of the 



physical parameters in (1.1 ). For smaller applied fields the energy of the minimizers per unit area 



in macroscopically large plates is 



Energy \ LAKH,^ f _ H^V f _ 
Area | A' + 47rM2 \ J \ ' ' ^ ' 

In particular, the energy per unit area of the plate scales as with the plate thickness and linearly 
(up to a slow logarithmic dependence) with the deviation of the applied field from the saturation 
field. This energy is achieved by trial functions consisting of periodic patterns of slender needle-like 
disconnected domains of magnetization opposing the applied field. In each unit cell of such a trial 
function the magnetization pattern refines toward the plate boundaries in a self-similar fashion (see 
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a) 





Figure 1: A sketch of the refining needle configuration, (a) The projection of one period of the 
domain pattern on the xia;2-plane. (b) A three-dimensional sketch of one period of the domain 
pattern. Shaded regions indicate the domains of magnetization opposing the applied field. 



Fig. [I]). This class of magnetization patterns is, therefore, a natural candidate for the precise form 
of the energy minimizers (see also Fig. [2]). 

In Section [4j we further investigate the asymptotic behavior of the energy in macroscopic samples. 
Under a physically reasonable assumption that the magnetization vector does not deviate strongly 
from the easy axis, we rigorously derive a reduced energy, whose minimum agrees asymptotically 
with the sharp interface version of the energy, see Theorem 4.1 The obtained result assigns a 



mathematical meaning to the /x* -method for computing the energy contributions away from the 
domain walls in a magnetization pattern, which was proposed more than half a century ago in 



the physics literature 33 . The obtained reduced energy, given by (4.2), practically coincides with 



that of an infinitely hard material in which the strength of magnetostatic interaction has been 
suitably rcnormalized. The latter explains why the behavior of minimal energy in both hard and 
soft materials is the same up to a certain factor in the macroscopic limit. Let us also note that 
for the same reason the energy per unit area becomes essentially independent of the saturation 
magnetization Mg in soft materials with fixed value of H^xt/Hs, see (1.2 1. 



In Sectionjsj we perform a formal asymptotic analysis of the reduced energy in (4.2 ) and establish a 
precise asymptotic behavior of the critical field Hc^ at which the only minimizer (global or local) is 
expected to be the uniform state, and the critical field Hc^ at which non-trivial minimizers emerge, 
see Theorem |5.1[ It turns out that asymptotically for macroscopically large plates 
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where Cq 



0.4368 and Ci 



Co, 



In 



(1.3) 



L^Mf [K + 47rM2 ) V AK 
0.5403. Ai H ^ Hc„ ^ the magnetization patterns are expected 
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Figure 2: Side views of magnetization domain patterns refining towards the boundary in bulk cobalt 
crystals: in the absence of the magnetic field (a) and in the applied field at 60% to saturation (b). 
From Ref. [f6]. 



to consist of slender, approximately radially-symmetric needle-like domains spanning the entire 
plate thickness and separated by large distances compared to the needle radius. Equation (f.3) is 
obtained from a reduced one-dimensional expression for the energy of needles, see (5.4 1. Solving 
the respective Euler-Lagrange equation exactly, we obtain the precise shape of the needle and, 
correspondingly, the expression in (1.3). 



Structure of the paper and notations 

The paper is structured as follows: In Section [2] we present the micromagnetic energy functional 
and introduce its sharp interface version. In Section [Sj we prove matching upper and lower bounds 
for bulk samples near the critical field. In Section [4j we derive a reduced model that captures 
the leading order energy in the macroscopic limit. In Section [5| we perform a further reduction of 
the energy and find the precise location of the transition to non-trivial minimizers by solving the 
reduced minimization problem exactly. 

We will denote a generic point in space by a; = (xi, X2, x^) = {xi, xa_), where xi is the component in 
the direction of the easy axis and x_\_ ~ (x2, ^3) is the component projection onto the plane normal 
to the easy axis. Similarly, we will denote the component of a vector v g in the direction of 
the easy axis by vi and its projection to the plane normal to the easy axis by v_l = (i'27i'3)- The 
spatial gradient is similarly separated into the components along and perpendicular to the easy 
axis: V = V_l). 

We use the symbols and > to indicate that an estimate holds up to a universal constant. For 
example A ^ B means that there are universal constants c, C > such that cA < B < CA. The 
symbols <C and 3> indicate that an estimate requires a small universal constant. For example, if 
we say that A < B ior e <^ 1, this is a short way of saying that A < CB holds for all e < £0 where 
Eo > is a small universal constant. By the symbol ~, we indicate asymptotic equivalence of two 
expressions: E.g. by writing A ~ i? for e ^ 1 and A ^ 1, we mean that for every S > 0, there are 
£0, Ao > such that \A/B ~ 1\ < S for all £ < £0 and all A < Aq. 

2 Physical model and sharp interface energy 

Micromagnetic energy: Up to an additive constant, the micromagnetic energy for a mono- crystalline 
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uniaxial ferromagnet (see e.g. 16 , using CGS units) is given by 



2./b3./b3 r-r' Stt ./5 



(2.1) 



Here, 17 C M"^ describes the region of space occupied by the ferromagnetic material, and the magne- 
tization vector M : ]R'^ — )■ satisfies |M| = Ms in Vl and M = outside. The terms in the energy, 
as they appear in the formula, are: 

1. The exchange energy favoring a uniform magnetization. 

2. The anisotropy energy favoring alignment of the magnetization with the easy axis. 

3. The Zeeman energy favoring alignment with the external field Hoxt- 

4. The stray field energy describing long-range Coulomb interactions of the "magnetic charges" 
V -M. 

5. A constant term, added for convenience. 



Also in (2.1 1, A is the exchange constant, Af^ is the saturation magnetization and K is the uniaxial 
anisotropy constant. The subscript "_L" denotes the components of a vector in the plane normal to 
the easy axis. 

Geometry of the sample: We consider a plate of constant thickness L, whose surfaces are oriented 
in the direction normal to the easy axis. To simplify the issues associated with the treatment of 
the lateral boundaries of the sample, we assume periodicity with period C in the plane normal to 
the easy axis. We hence write il = (0,L) x T, where T = [0,£)^ is a torus with periodicity £. 
The periodicity assumption, however, is not essential, as long as the energy of the minimizers is 
extensive in £, i.e. we have inf £ = 0{C^) as £ — >■ 00. As we will show below, this will indeed be 
the case. Also, the external field is assumed to be in the direction of the easy axis (and hence also 
normal to the material surface): 

Hext = ^^oxt ei , 

where ei is the unit vector in the direction of the easy axis. 

Rescaling: As usual, we first introduce the exchange length Zcx = ■\/j4/(47rA/|) and the dimensionless 
quality factor Q = K/{A-kM'^) [16]. Introducing m = M/M^, hext = Hcxt/(47rM5), £ = £/L, and 
measuring lengths and energy in units of L and , respectively, we can then rewrite (2.1 ) as 

£[xn] = [ |||Vmp-fQ / [m^p - / (2h,xt • m - h • m , (2.2) 
Jn Jn Jn 

where 17 = (0, 1) x T is the rescahng of fi, with T — [0,£)^, and |m| = xn, where xn denotes the 
characteristic function of O. The dimensionless stray field h is defined as the unique (see e.g. [6]) 
solution in L'^{R x T; R^) of 

V X h = 0, V • h = -V • m in M x T, (2.3) 
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where (2.3) is understood in the distributional sense. 



Sharp interface energy: In a bulk uniaxial material, the magnetization is expected to lie mostly 
in the direction of the easy axis, i.e. m ±ei. The regions with different orientations of the 
magnetization are expected to be separated by thin Bloch walls 16 p. 215]. A Bloch wall is 
characterized by a transition layer of thickness w « Q~ s 1^^ in which the magnetization rotates in 
the wall plane, thus avoiding the creation of a stray field. Taking advantage of the observation 
in fl9| p. 367], one can directly estimate the anisotropy and exchange terms in the energy from 
below for any (5 > fixed as 



/ ( 



^|Vm|2 + 0|m^|2 



> £ 



|Vmi| 



(2.4) 



with the notation m = miei + m_L and where we have introduced 



L 



(2.5) 



In fact, the one-dimensional Bloch wall profile attains equality in (2.4) (see |16|), which implies 



that the term in the right-hand side of (2.4) should actually well approximate the term in the left- 



hand side for the energy minimizers. The condition of validity of this approximation is that the wall 
thickness w remains much smaller than the characteristic length scale of the magnetization pattern. 
In particular, one should have w L. This condition is achieved for sufficiently thick plates. In 
fact, large thickness is also a necessary condition for branched domain patterns to be observed in 
ferromagnetic materials 16 . Therefore, in the present context one is naturally interested in the 
asymptotic behavior of energy for large values of L or, equivalently, in the limit e — !■ with all 
other dimensionless parameters fixed. 



Dropping the gradient term in (2.2) where |mj^| < S and combining it with (2.4), one can see that 



£[m]> e [ iVmil+g/ \m^\^ ^ [ {2h 

J\in±\>S J\m±\<S Jn 



7^ ^ 



(2.6) 



This motivates the introduction of a sharp interface energy, in which the gradient-squared term in 



(2.2) is replaced by the total variation of mi (see also [6j[7|). We note, however, that the sharp 
interface energy is basically a tool to approximate the behavior of the full physical energy in (2.2) 



and, therefore, can be tailored to our advantage. We choose the sharp interface energy in the form 



\ra±\^+S^Q \ni±\^- (2hext • m + h • m - /z^^^) . 

(2.7) 



Here < (5 ^ 1 is an arbitrary "cutoff" parameter, whose precise value is inessential (hence the 
index 5 is dropped from the definition of E), and is the truncated version of mi: 



1-^2^ mi{x)>l-6^, 
mi{x), -1 + (5^ < mi(.T) < 1 - (5^^ 
-1 + S^, mi(x) < -1 + 5^. 



(2.8) 
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The advantage of using mf in (2.7) instead of mi is that, consistently with (2.6 1, the interfacial 
term does not contribute to the energy away from the domain walls, where mi « ±1. Importantly, 
E provides an ansatz-free lower bound for £: 

£[m]>{l-5^)E[m\, (2.9) 



which can be easily seen by retracing the arguments leading to ( 2.6 1. Furthermore, since in the limit 



£ — >■ the transition regions between different directions of m are expected to become 0{eQ ^) 



thin and the inequality in (2.6) to become an equality for minimizers, in view of arbitrariness of 5 
one should expect that inf£ ~ inf _E for e ^ 1. In the following, we will prove that this relation 
holds in the sense of scaling, i.e., for sufficiently small e > 0, we have mi£ ^ inf E. 

Critical external fields: Clearly, when the applied field h^xt is sufficiently large, the minimal energy 
configuration will be such that all magnetic moments are aligned with the field, i.e. m = ei^n- 
For smaller external fields, the niinimizer is attained by other configurations. The external field 
strength, at which the uniform magnetization m — eixn looses its optimality is denoted by hc^ . Let 
us also note that appearance and disappearance of patterns as a function of the control parameter 
in systems of this kind is often accompanied by a hysteresis. Therefore, non-trivial critical points 
of the energy may persist even for fields larger than /i^ . The critical field at which these critical 
points disappear will be denoted by hc^. 

The saturation field hs is defined similarly in terms of the relaxed energy. In our setting, it is the 
variant of £ where the surface energy is not penalized, i.e. 

£rci[m]= / Q|m_Lp- / (2hext • m - h • m + /i^^J . 
Jn Jn 

The set of admissible functions for £^.^1 is given by all m satisfying |m| < 1 in fi. This relaxed 
constraint in the above calculation can be justified by looking at small-scale oscillations of m. It 



is related to the fact that frci is non-convex, see e.g. 10 21 . One expects that the relaxed energy 



gives the leading order behavior of the minimal energy, i.e. 

lim inf £\m\ ~ inf £ioi[ml- 

£^0|m| = l |m|<l 

The saturation field strength hg is defined as the field strength at which m = eixo looses its 
energetic optimality in terms of the relaxed energy £^ci and is expected to be close to when 
£ < 1. 

To understand better the behavior of £roi, let us first introduce the notation for the average of a 
quantity / — f{xi, X2, X3) over T at fixed xi e M. We use the notation f{xi) :~ -p Jj, f{xi, •). We 
note that the solutions of ( |2.3[ ) have the following basic properties (the proof is by an elementary 
integration by parts): 



Lemma 2.1. Let h e L^{R x T;E-^) be a solution of {2.3). Then 



hi{xi) = —mi{xi) for a.e. Xi G M, (2-10) 
hm= f |h|2, / |hp < / Imp. (2.11) 
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Using (2.10) and (2.111, one easily computes that 

i |2.ioy2.ii[ 



£rcl[m] 



(/icxt + hi) 



where we have used Jensen's inequahty in the second inequahty. Hence, 



frci[m] > inf / (/lext - mif = \ 

|m|<Uj7 (1 



< /loxt < 1, 



^cxt)^, ^cxt > 1- 



(2.12) 



(2.13) 



On the other hand, equahty in ( 2.13 ) is achieved by using the trial function m = min{l, h^xt} ^iXn- 
Therefore, m = eixn is the minimizer of (2.131 if and only if /icxt > /is = 1, which is precisely the 
saturation field. 

Since we are interested in the bifurcation from the uniform to a patterned magnetization occurring 
near saturation, we introduce a parameter A which measures the deviation from saturation: 



1- A, 



(2.14) 



where < A <C 1 means the applied field is just below the saturation threshold. One question we 
want to address is how to calculate Ac^ and A^ corresponding to the critical fields hcg and /ici ■ 

Reformulation of the sharp interface energy: We now derive an expression for energies £ and E in 
new variables which make our analysis more convenient. We introduce 



u = m- (1 - A)xoei, 



V = h + (1 - A)xoei. 



(2.15) 



Then, using Lemma 2.1, ( 2.14 ), ( 2.11 ) and ( 2.15 ), one gets 

/(/i,\t-h.m-2hexfm)= / |v|2 - 2(1 - A) / ( 
Jn Jrxt Jn 



Vi + Ui) 

Therefore, we can rewrite the energy E from (2.7) as follows (with a slight abuse of notation, we 

|u±P+/ |vP + <52q/ |u^|2, (2.16) 



view E from now on as a function of u instead of m) 



E[u]=e I iVufi +Q 



/{|ui|<5} JKxT J{|ui|>5} 

where, as before, ui and uj^ denote the components of u along and normal to the easy axis, 
respectively, uf ~ m\ — 1 + A, and v solves 



-V(/3, A(/7 V • u 



X T, 



(2.17) 



The set of admissible functions for (2.16) is given by 
= {u e BV{M. X T;M3) 



u+ (1 - A)xf2ei| = xn}- 



Similarly, the expression in \2.2\ can be rewritten as 

-2 



4Q 



iVul- 



Q 



(2.18) 



For simplicity of notation, we take the same admissible class A for £ as well, setting £[vl\ — +oo, 
whenever U|j7 ^ H^{n). 
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3 Scaling of the energy in bulk samples 



In this section, we investigate the scaUng behavior of the energy of minimizers in the case of bulk 
samples corresponding to the limit e — > 0. The main part of this section will be concerned with the 



sharp interface energy E defined in (2.161. The connection to the diffuse interface energy £ is then 
shown in Section 13.51 



The model has three dimensionless parameters: e, A, Q. In particular, we are interested in the 
case of macroscopically large samples near critical fields, i.e. e ^ 1 and A ^ 1. Our result shows 
that for sufficiently small e and A with fixed Q there are exactly two different scaling regimes, 
each corresponding to a particular pattern of magnetization attaining the minimal energy scale. 
Introducing 

we have the following result for the sharp interface energy E: 

Theorem 3.1. Let A < In Ap and ^>7^3e5A~5|lne|^3. Then for e <C 1 and A <C 1, we have 



inf E\u] ^ min/A^, 73£3A|lnA|3V 



The first regime corresponds to a uniform magnetization along the applied field, while the second 
regime is achieved by branched magnetization patterns. Note that as long as Q > 1, the particular 
value of Q does not affect the scaling of the minimal energy. This indicates that the restricted 
model corresponding to Q — oo, i.e. when m = ±eixn, captures the essential features of the 



general model in (2.16). On the other hand, for Q 1 the effect of anisotropy only has the effect 
of renormalizing the minimal energy scaling by a factor of Q ^ . This will be further discussed in 
Sec. |4]with the help of a reduced sharp interface model. 



Combining the results in Theorem 3.1 with (2.9) and the constructions of Sec. 3.5 the full micro 



magnetic energy £ satisfies the same scaling 

Theorem 3.2. Let A < 7^| In Ap and £ > -f~i ei X^i\\n£\^i . Then for £ <C 1 and A <C 1, we have 

— inf £[u] - -^"^ ^ ^ 75£5|ln£|5, 

P ueA |73£3A|lnA|3 for A > 73£3|ln£|3. 

This theorem implies that for small enough values of A(£) the minimal energy scaling is achieved 
by uniform magnetization pattern (the monodomain state: m = Bixn), while for sufficiently large 
values of A(£) the optimal energy scaling is achieved by a branched domain pattern, as £ — > 0. The 

1 2 I I 1_ 

transition occurs at Ac ~ 73£3 1 ln£| 3 . 

The analysis techniques we employ in this section go back to the work of Choksi and Kohn in 
[6][7], who analyzed the energy of ferromagnetic plates in the absence of a magnetic field. In our 
analysis we identify the optimal dependence of the minimal energy on the parameter A which is 
not addressed in (6][7]. In our analysis, we also apply tools from related works in the framework 
of type-I superconductors [5^. There the authors derive the scaling of the energy for the type-I 
superconductor near critical field. We note that the super conductor model is more rigid, since there 
the two different phases are described by the characteristic function x which only takes the discrete 
values and 1 and a divergence free magnetic field i?, whereas in our model the magnetization m 
is allowed to take all values on the unit sphere. 
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3.1 Preliminaries 



In this section we collect some useful results before addressing the proof of the upper and lower 
bound in the next two sections. 



Control on (2.16) yields information about ui and vi on each slice. As expected, the stray field 
favors zero average of ui on each tangential slice: 

Lemma 3.3. Let u e A and let v e L^(M x T;M^) satisfy ( |2.17[ ). Then for every Ci > and 
C2 > there exists a constant c > 0, such that if E[\\\/f' < cA^ th 



en we nave 



Jo 



Furthermore there exists I C (0, 1) with |/| > 1 — C2, such that for all a ^ I 

\ui{a)\ < ciA. 



(3.2) 



(3.3) 



Proof. We first note that in terms of u and v, in view of (2.10), we have vi — —ui for a.e. a S 
By Jensen's inequality, it then follows that 

['\uifdx, - / \u,f = [ |v|2 < / |v|2 < E[u], 
Jo Jn Jn Jn 



and (3.2 1 follows. Inequality (3.3) follows from (3.2) by an application of Fubini's Theorem. □ 

The main ingredient for the proof of the lower bound is an estimate that characterizes the transition 
energy, i.e. the cost for the magnetization to vary between a tangential slice {a} x T and its value 



zero outside of the sample. The idea to estimate such transition energies was introduced in 20 
and has been subsequently applied also in e.g. [sllsl. 



Lemma 3.4 (Transition energy). For every u € ^ there exists I C (0,1) with \I\ > |, such that 
for all a G I and for all if) £ H^{T), we have 



{a}xT 



2[u] (7 5||Vi/'||l2(x) + UWmT) 



Proof. Let us first assume that u S Cj?°(M x T), with u = outside of [— 1, 2] x T. Let v be defined 
by (2.3 1. Noting that ip does not depend on xi and using integration by parts, for any a G (0, 1) 
and b £ (—2, —1), we then get 



b JT 



di{ui{xi, dxi 



pci p pci p pa, p 

/ "^^-u^ixi,-) ipdxi- Vj_-vj_{xi,-) ipdxi- 

Jb JT Jb JT Jb JT 

f f u^{xi,-)-W±tl;dxi+ f f 

Jb Jt Jb Jt 



diVi{xi, •) ipdxi 



v_l(2^i, •) ■ c?a:^i + / vi{a, ■)ip dxi — j vi{b,-)^. 

(3.4) 



11 



By Fubini's theorem, there exists b E (—2, —1) and / C (0, 1) with |/| > ^ such that for aU a E I, 



l«il 



{a}xT 



2 < 



{6}xT 



(3.5) 



inequahty and (2.16). 



The statement then follows for all a G / from (3.4), (3.5) and by application of Cauchy-Schwarz 



Now consider a general u G In this case, u can be approximated by a sequence of functions 
E C^(R X T) such that = outside of [- 1, 2] x T and such that u-^' ^ u in L2(k x T;R^) 
/rxtI^uI' see [mI 



xT 



iVu^ 



By (|2.1lD, we also have v in L'^{R x T; M^), where v-'' 



and /k 

denotes the stray field of . Taking a subsequence, if necessary, we also have convergence u-' — !• u, 
v-' — )■ V in L^{{a} x T;M'^) for a. e. a G M. Using this approximation, the lemma follows. □ 

We will also use the following technical lemma of De Giorgi (see, e.g., [5| Lemma 3.1]): 

Lemma 3.5. Let S CZ T be a set of finite perimeter, and let r > be such that r\dS\ < ^15*1. Then 
there exists an open set S d T with the properties 

(i) There is a considerable overlap of S with S, in the sense 0/15*0 5*1 > ^IS"!- 
(li) For all t > 0, the set := {p E T : dist {p,S) < t} satisfies \S*\ < |5|(1 + (t/r)^). 

3.2 Ansatz-free lower bound 

In this section, we present the proof for the lower bound. We need to show: 
Proposition 3.6 (Lower bound). For e ^ 1 and A <C 7, we have 

inf E[u] > minlA^, 75£3A|lnA|H. 

Proof. Following the ideas in [s], we argue as follows. Recall that in view of Lemma 3.4 the energy 
is bounded below by a Sobolev-type norm of negative order on u, evaluated on a generic tangential 
slice {a} x T. In this proof, we combine this with control that we have on the surface energy and 
anisotropy energy on a generic slice. The proof is divided into five steps. 

Step 1: Identification of tangential slice. We will argue by contradiction. Hence, we may assume 
that the energy does not satisfy the lower bound, i.e., there exists u E A, such that 



E[\i] < ^^minjA^ 75e5A|lnA|5| 



(3.6) 



for some e ^ 1 and A ^ 7. We choose a E (0, 1) such that the assertions of Lemma 3.3 and Lemma 



|3.4| hold. By Lemma |3.3| and by Fubini's Theorem, we may then assume that a is furthermore 
chosen such that 



|ui(a)| < A 



IVuil 



{a}xT 



^ — t mm ■ 



73 A| In A| 3 



£3 



< 



E H 



£ min ■ 



{a}xT 



A^ 73e3A|lnA|3 



(3.7) 
(3.8) 

(3.9) 
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Step 2: Structure of magnetization. We next analyze the magnetization on the shoe A :— {a} x T in 
more detail. In view of the upper constructions, we expect that the regions where the magnetization 
points in the negative xi-direction are small needle-shaped domains. The restriction of mi to the 
slice, therefore, is expected to be negative on a number of small circular domains. In the following, 
we give a precise version of this heuristic picture. We define the set (where m points "to the 
right", i.e. in direction of ei. The notation "to the right" is in accordance with the figures) and 
the set A- (where m points to the left) by 



A+ := {xeA: < ui < A}, 



A^ := [x e A : -2 + A < ui < -1 + A}. 



It is also convenient to define a transition region Aq (where m also points to the right) 

Ao:= {xe A : -l + \<ui< ^A}. 
We first note that when mi > 0, we have \u±\'^ = \m±\'^ = 1 — ml > 1 — mi = A — wi, i.e. 



|u_l|"' > A-iii 



in Aq U A^ 



(3.10) 



We claim that the region A^ of "reversed magnetization" is concentrated on a small set with total 
area of order X£'^, and that the transition region is even smaller. More precisely, we claim that 



1^^ 



\Ao\ « 



\A. 



(3.11) 



Indeed, by (3.101, we have |uj^| > 4 in Aq, and so 



< 



,2 < A ^2 

' AQ ^ ' • 



Here we also have used that by assumption Q > 7 ^ A. Choosing ci = in Lemma 



view of I A 



01 



\A 



+ 1 



-^A^2 

16 



ui < {-l + X)\A^\ + X{\Ao\ + \A+\) = -\A^\ + X£^ 



r - |A_|, we get 

^ r 

~ J A 

SO |A_| < i|A^2 ^ £2 g^j^^j^ therefore, \A+\ > ^i'^. Similarly 



3.3 



and 



16 



Ar 



- f ui < 2\A_\ + f iX-ui)-'^X\A+\ 

J A JAo 

2|A_|+ / \n^\'~lXf< 2|A_ 



CE 



\Xf 



2\AA~\X^\ 



Hence > ^XP , and so Xt"^ . This concludes the proof of (3.11). 

Step 3: Identification of a regularized region. In the previous step, we showed that the reversed 
magnetization region A_ occupies a small fraction of A. By the upper bound constructions of 
Lemma 3.8 below, we would expect that A_ is divided into a controlled number of similar size 
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circular domains. The construction also suggests that in the core region of the plate, the typical 
radius r of these circular domains and the typical distance a between them are given by 



£3 



£3 



73 1 In A| 3 



7 3 A 2 1 In A 1 3 



(3.12) 



In the following, we use the co-area formula and the isoperimetric inequality to get a rigorous 
variant of the above heuristics. We replace A- by a larger set S with A- CSC A- U Aq. The 
reason to choose S instead of A_ is that we cannot exclude a concentration of surface energy on 
OA-. We claim that there is c G (^57 ^i) such that the set 



satisfies 



\dS\ 



S := { X e A : m < c} 

\s\ ^ \e, 

2 . /a2 7^A|lnA|i 
I mm < — , 1 

e £3 



(3.13) 

(3.14) 
(3.15) 



Indeed, by ( |3.13[ ) it follows that A_ C 5 C A_ U Aq and ( |3.14[ ) follows by ( |3.1l| ). Furthermore, by 
the co-area formula 

/ iV^uil = / n^ix^ e A : wi(.T^) = t}) dt, 

J A Jr 

and Fubini's Theorem, there exists c G (— |, — |), such that S satisfies 



\dS\ 



< 



|V_LUi|9 e'^ mm- 



A2 73A|lnA|s 

£ £3 



Estimates ( |3.14[ ) and ( |3.15| together with ( |3.12[ ) yield 

r \dS\ < \S\. 



(3.16) 



Heuristically, the estimate in (3.16) means that S, roughly speaking, splits into a collection of disks 
of diameter much larger than r. This disagrees with the expectation from the upper construction 
and will lead to a contradiction. 

We next replace S by another set S, still satisfying all the relevant properties of S. Additionally, it 



grows in a controlled way upon "thickening". More precisely, in view of (3.16) and by Lemma 3.5 
there is a set S with 



\sns\ > ^\s\ 



and such that for alH > r, the t-neighborhood S* of satisfies 



(3.17) 



(3.18) 
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Step 4- Definition of a suitable test function. We now define a logarithmic cut-off ip £ H^(T), witli 
< V' < 1, around S. Let 



il^ix^) := (/^(dist S")), 



wfiere 



^(t) := 



1 for < i < r, 

forr<t<a, 

m(a/r) — — ' 

for a<t. 



(3.19) 



with r and a defined in (3.12). A direct computation, following then yields 



< 



|lnA| 



and 



< 



i2|lnAr 



(3.20) 



For the reader's convenience, we show the first estimate in (3.20), the proof of the second inequality 
proceeds similarly. Since tp = 1 on S, and since \dS*\ — -^\S |, we have 



^ ^ \S\+ ip{t) \dS*\ dt 



^'(t) \S'\dt 



xe 



''2 1 In a/r 



t dt < 



r-2|lnA| 



llnAI 



Step 5: Proof of the lower hound. We are ready to give the proof of the lower bound. It is based 
on application of Lemma |3.4| and on a duality argument, using the test function -0. We claim that 



Xf < 



Uiljj. 



(3.21) 



Since ■0 = 1 and ui < — ^ in S', and since > and ui < in it follows for A ^ 1 that 



uiip 



Ui0 



sns 



A\{SnS) 



1*10 > -|S'nS'|-A / V 



l |3.14tj3Tfl l 



xt 



Application of Lemma |3.4| then yields 



Xf 



We hence obtain 



Mi0 < S3[u] I ^||VV'|U2(t) + ||0I1l2(t) 

V72 



£;5[ui 



70^1 In A I 



E[u] > fmm{jX^a^\lnX\,X^} ^ ^Vrin {75ei As | In A|, A^}, 



contradicting (3.6). This concludes the proof of the Proposition. 



□ 



3.3 Sharp interface constructions 



In this section, we present constructions that achieve the optimal scaling in Theorem 3.1 We have 
two different regimes. For smaller values of A, the optimal scaling of the energy is achieved by a 
uniform configuration, while for larger values of A, the optimal scaling is achieved by a self-similar 
structure. Note that the constructions in ^ do not yield the optimal energy in the case of near 



15 



saturation field. Instead, our constructions are an adaptation of constructions introduced in [s] for 
a model of type-I superconductors. While our constructions have a similar self-repeating structure 
to that of [g], the definition of the involved functions is different due to the constraint |m| = xn in 
our model. Our constructions are also different from those of [Sj in the geometry of the magnetic 
domains, and are intended to better mimic the behavior of the minimizers. The main result of this 
section is: 

Proposition 3.7 (Upper bound). Suppose that A < 7^|lnAp. Then for e <C 1 and A <C 1, the 

scaling of the minimal energy E is bounded above by 



I; inf Ehi] < min(A^, 73e3A|lnA|3l 



£2 

Let us remark that the logarithm in the scaling of the energy is a consequence of the fact that 
the leading order contribution of the stray field energy is given by interaction on tangential slices, 
where the stray field potential behaves logarithmically. 

We first note that by choosing the uniform magnetization u — Aei^n, we immediately recover 
the upper bound ^ infugy!i i?[u] < A^. The cross-over to the branched regime occurs at A ^ 
j3£3 1 In A| 3 . It hence remains to construct an optimal upper bound, if A is larger than this thresh- 
old. In the remaining part of this section, we present such a construction. The corresponding 
estimates are then given in Section |3.4[ thus completing the proof of Proposition |3.7[ 



Before going into the details of our constructions, however, let us recall some other constructions 
that have been proposed in the literature over the years (for simplicity, we will only discuss the 
case Q ~ 1). The first estimates of the minimal energy for the bulk uniaxial ferromagnets go back 



to the work of Landau and Lifshitz 24 and Kittel [18]. Those constructions were proposed for 
zero applied field. In fact, the Landau-Lifshitz construction cannot be easily extended to the case 
when the domains opposing the applied field occupy only a small volume fraction of the sample. 
The Kittel construction, on the other hand, can be modified to account for small volume fraction. 



resulting in an energy scaling Ei~ ^ e^\\\nX\ 22 . However, since it consists of a "striped" 
domain pattern, its interfacial energy turns out to be too high at small A. This issue can be 
addressed by modifying the geometry of the domains into a lattice of cylindrical "bubbles" , whose 
energy my be estimated as E£~'^ ~ e^A [s 13 . A comparison of this estimate with the result of 



Proposition 3.6 shows that, although the bubble construction provides a slight improvement over 
the stripe construction, it is highly non-optimal in its scaling behavior with respect to e. We note 
that, in fact, any domain configuration, in which the domain walls are aligned with the easy axis 
cannot do better in terms of energy, and so branching is inevitable for sufficiently small e to reduce 
energy On the other hand, if a tree- like branched domain structure is used (see [sj Sec. 4.3] for 
the construction in the case of type-I superconductors), it is not difficult to show th at th e energy 
will scale as Ei~'^ es As . Once again, comparing this with the result of Proposition 



3.6 



one sees 



that, while the considered configuration gives the optimal scaling in terms of the dependence of the 
energy on e, it is highly non-optimal in terms of A. These observations indicate that the minimizers 
of E may not have the geometric characteristics of any of the domain patterns considered above 
when e <C 1 and A ^ 1. In the following we present a construction which achieves the scaling in 



Proposition 3.6 thus demonstrating that this scaling is optimal. 



We begin by fixing the basic geometry. The geometry is an adaption of a recent self-similar con- 
struction for the type-I superconductor model |5j. However, contrary to the construction in 5 , 
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closure transition cells core cells transitiorr cells closure 

cells cells 



Figure 3: A side view (in the a;ia:2-plane) of a sample partition containing 16 core cells with 4 
generations of refining cells on each side. 



our construction includes closure domains. These closure domains are, in particular, necessary to 
achieve the optimal scaling of the energy in the case of soft materials, i.e. Q 1. Based on this 
geometry, we construct two different magnetization configurations u"*^ and u"^^. The first config- 
uration u'^^ avoids anisotropy energy entirely and is optimal for Q ^ 1- The second configuration 
u'^-'^ avoids most of the stray fields and is optimal for Q < 1. 

Sample geometry. We divide the material plate into three spatial regions: the core region, 
the transition region and the surface region, see Fig. [3l These regions form 5 layers symmetrically 
with respect to the plate's mid-plane. By this symmetry, it is, therefore, sufficient to describe the 
constructions only in the left half of the sample, i.e. ior < xi < ^. 

The core region is partitioned into equal rectangular cells with height ho in the normal direction 
(xi-direction) and length oq in both tangential directions. These cells are adjacent on the left to 
a system of M layers of self-similar cells in the transition and surface regions that refine from the 
core region towards the boundary (see Fig. [3| . Each generation of cells is described by its height 
hj in the normal direction and its extension aj in both tangential directions, with j — 1, . . . , M. 
We also define a parameter rj — aj /\/2Tr, which will be the maximum needle radius in the j-th 
generation of cells. In every generation, the width of the cells decreases by a factor 3, i.e. 

and rj+i = ^. (3.22) 

In particular, the number of cells is multiplied by a factor 9 in each new generation. The algorithm 
is terminated after M iterations. We will specify oi, {hj}, and M in the sequel. In particular, these 
parameters will be chosen, such that the union of all cells exactly covers il, i.e. 

M 

ho + 2^hj = 1. (3.23) 

We differentiate between core cells (which for simplicity of notation we identify with generation 0), 
transition cells (generations 1, . . . , M — 1) and closure cells (generation Af ). 
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Figure 4: Refinement towards the sample boundary (two generations are shown). 

Globally, the geometry of the construction consists of a collection of needles rescaled to fit into the 
collection of cells just constructed, refining in the direction of the boundary, see Fig. [4] We capture 
the region occupied by the needles by the characteristic function x G BV{9? , {0, 1}), which will be 
specified in the sequel. We thus give the definition of % on a rescaled cell 

Z = [Q,h]y.K, ^=[-f'f]^ (3-24) 

with height h and width a. We furthermore denote the tangential boundary of the cell by d^Z := 
[0, h] X dK. The corresponding "maximum needle radius" r is defined by 




(3.25) 



The definition of x o^i ^-^ly cell with arbitrary extension (in the left side of the sample) is then given 
by a rescaling of this cell. 

Geometry of a transition cell. Consider a transition cell Ztms first. The cell geometry is 
characterized by nine needles, see Fig. [sjja). The largest needle is located in the center of the cell 
and grows into positive xi-direction, while the other needles are smaller and grow in the negative 
xi-direction. All needles are axially-symmetric around their corresponding center lines given by 
— a;||^\ i = 1, . . . , 9. The large needle is located in the center of the cell, i.e. = 0. The radii 
of the needle cross-sections on tangential slices are functions of Xi. The radius of the large needle 
is denoted by pi{xi) := p+{xi). The radii of the 8 small needles are given by pi{xi) := p_(xi) for 
i = 2, . . . 9. The characteristic function x is defined by 

9 

x(xi,a:_L) := ^ i7(pj(a;i) - |a;_L - |) for (xi, a;_L) G ^'tms, (3.26) 

i=l 

where H is the Heaviside function, i.e., H{s) = 1 for s > and H{s) = for s < 0. It remains to 
specify the radii p± for the large and small needles. 
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Figure 5: Geometry of a unit cell: a) transition cell; b) closure cell. 



At the tangential faces, i.e. at xi = and xi = h, of the cell the needle radii are defined by 

p+(0) = p_(0)-ir, P-{h)^0, p+{h)^T, 

respectively. This means that at a;i = all needles have the same radius. At xi = h, the radii of 
the small needles are zero. We choose p±, such that throughout the cell, the cross-sectional area of 
the needles is constant, i.e. 



Trp1{xi) + 87rp^(a;i) 



(3.27) 



To avoid further complicating the constructions, we assume that the profile of the small needles is 
conical at the tip: 

hp_{xi) = r{h-xi), for (1 - a)/i < xi < /i, (3.28) 

where < a ^ 1. The precise value of a is inconsequential and, in fact, for the zero anisotropy 
configuration u"^^ could even be taken to be zero. In view of (3.27), equation (3.28) defines p+{xi) 
for (1 — a)h < xi < h. Furthermore, on most of the cell, i.e. for < xi < (1 — a)h, we choose p+ 
to be the linear interpolation connecting the values of at = and xi = (l — a)h. In turn, 
P- is defined by (3.27). Note that as a consequence of the conical profile at the tip of the small 
needles (see (3.28)), it follows that 



\p'4x,)\, |p'_(xi)| < ^ 



for all xi e [0, h]. 



(3.29) 



Geometry of a closure cell. In any closure domain cell Z^s, there is only one large needle 
along the center of the cell, seejsjb). The radius /5cis(a;i) satisfies the conditions Pcis(O) = and 
Pcis{h) = r. Analogously to the above construction for the transition cell, we choose a needle shape 
with a conical tip. Define 

x{xi,xj_) := H{pcisixi) - \x±\) for x e Zeis. (3.30) 

In particular, we also have 

|Pcis(2^i)l < T foraUxi G [0,h]. (3.31) 
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Geometry of a core cell. In each core cell ^corc the function x is assumed to be a characteristic 
function of a straight cylinder with the radius equal to that of the needle on the adjacent side of the 
transition cell. The overall geometry of the magnetization pattern for one core cell and 3 refining 
generations is presented in Fig. [T] In the sequel, we give two different magnetization configurations, 
u'^^ and u"^^, based on the geometry described above. 

The magnetization u'^-'^. We first define the anisotropy-free configuration u"^-*^ by 



AF 



{X-2x{x))ei 



ill Zqqj-q^ ■^trnsi ^cls-i 



(3.32) 



and zero outside SI. Note that the stray field of u'^^ is created by surface charges on the needle 
interfaces and at the sample surface. We next define an auxiliary field v"^^. In every transition cell 
■^trns, we define 



where ip is a solution of 



in Zti-r 



and 



in Ztrr 



du^Lp = on d±Ztrns- 



Note that v"^^ is uniquely defined, since by (3.261 and (3.27) we have 
uf^^ (xi, xj_) dx 







Va;i e [0,/i], 



(3.33) 
(3.34) 

(3.35) 



K 



implying the solvability condition for (3.34), in view of the fact that u"*^ = u^^'^ei. Also note that 
the corresponding field v'^-'^ is the approximation of the stray field assuming that magnetostatic 
interactions in tangential slices are dominant. 



In closure domain cells, (3.351 does not hold any more. Hence, in this case we define v'^^ by 

hi Zeis, 



~AF 



= -Vj_tpi - diip2 ei 



(3.36) 



where W _]_^pi approximates stray field interaction in tangential directions, while diip2 approximates 
stray field interaction in the normal direction. We define (pi as a solution of 



A_L(y5i = V-u ^ ~ ^'j^V-\i^^{xi,x±^dx±^ in Zeis and d^^'Pi 
The function ip2 is defined as a solution of 



on 9_L^cis- (3.37) 



dfip2{xi) = Ik^ ' xj^) dx± in Zds and diip2 = on {xi = 0} x K. (3.38) 

Similarly to ( 3.34[ ), up to a constant, problems (3.37) and (3.38) are indeed uniquely solvable for 
(pi and ip2- Finally, we set v"^^ = in the core cells Zcorc and outside J7. 

The magnetization u'^^. Here we construct an approximately stray field-free magnetization 
u'^-'^. Although it would be natural to consider u — u"^^ + v"^^ as a trial function, in view of the 
fact that V • u = in this case, this function is not admissible, i.e. it does not belong to A. For 
this reason, using the construction of u'^^ and above, we define an auxiliary function v"^^: 



~SF 



{x):= (l - ^1 - |v^^(x)P )(u^^(x) + (1 - X)xnix)e,) in Zt^s, 
ix):= (l - '^1 - l^t^ixW ) (u-^^(^) + (1 - A)xo(x)ei) + v^''ix)e, in Z,i,. 



(3.39) 
(3.40) 
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As will be shown in Sec. 3.4 this definition is well-posed, since in our construction j-v"^^! <^ 1. We, 
therefore, set 



,SF 



(3.41) 



It can be easily checked that by our definition of v"^^, we have u"^^ e A. Furthermore, 
constructed to have small stray field (see Sec. 3.4). 

Localization of the stray field. Note that our constructions are such that 



is 



V • u^^ + V • = 



and 



V • u 



SF 



V • v^^ = 0, 



(3.42) 



in M X T. We will use this information to localize the estimates for the stray field energy (for the 
original idea in the context of ferromagnets, see [6]). Let us note that for every vector field u and 
its stray field v (in the sense of (2.17)), we have 



Ivp = inf 



(3.43) 



where the infimum is taken over all fields v € L^(]R'^; M'^) satisfying 

V-v + V-u = (3.44) 
distributionally. This motivates to define for any u ^ A and for any v e L'^(M x T; M^), the energy 



^[U,V] =£ / \Wui\+Q / |Ui 



(3.45) 



Hence we have E[u] < £'[u,v], whenever v satisfies (3.44 1. In view of (3.42) it then follows that 



(3.46) 



The advantage of the quantity E is that it is local in both of its parameters. We hence define the 
restriction of E on any set A C M x T by 



^l^[u,v] =e / \Vui\ + Q / |uj. 



Ann 



Ann 



3.4 Estimates 

In this section, we give the estimates corresponding to the branched needle construction described 
in the previous section, thus completing the proof of Proposition |3.7[ We will show that 

Lemma 3.8 (Needles). Forj^es |lne|s < A < T^llnAp, we have 

1 



inf £^[u] < 73e3A|lnA|: 



(3.47) 
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In view of (3.46), it is enough to give the estimate (3.471 in terms of the locahzed energy E, defined 
in (3.45), using the approx imat e stray fields v"^^ and v"^^ defined in the previous section. Before 
giving the proof of Lemma 3.8 we estimate the restriction of E onto a single transition or closure 
domain cell. We note that as long as Hq is bounded away from 1, the estimates for the core cells 
trivially result in the same upper bounds as for the transition cells. Therefore, in the following we 
do not include explicit arguments for the core cells. 

We first give the estimates for u'^-'^: 

Lemma 3.9 (Energy of a transition cell for u'^^). Consider a transition cell .^trns with height h 
and maximum needle radius r. Suppose that the needle is slender, in the sense of r <^h. Then 



E, 



K^v^"^] < erh 



InAI 



where u is defined in (3.32) and v is given by (3.33). 

Proof. By the definition of u"*^ , we immediately get the following estimate for the surface energy, 

S,AF\ 



: I IVw^'^^l ^ / |Vx| < erh, 



where we used the assumption r <^ h. It remains to give the estimate for the stray field part of the 
energy. In view of (3.32), we get 



V • u 



AF 



dmt^ = -2Y,5{p,{xi)-\x^-xf\)d^p,{x^), 



(3.48) 



where 5 is the Dirac (5-function. We expect a logarithmic blow-up near each of the needles. This 
motivates to decompose defined in (3.33), by 1^9 = (^^'^^ + ip'^^\ where (/j^"^ is given by 



ip^°\xi,Xi_) = ^2pi(a;i) aipj(xi)iJ(|a;^'' - - pi(a;i)) In- 



(3.49) 



»=i \x^-x^l>\ 

and where, as before, H is the Heaviside function. Let us for the moment assume that the leading 
order contribution to the stray field energy is due to (^(°). In view of (3.49), it is easy to calculate 

r^llnAI 



and 



r 
h?' 



(3.50) 



We then get 



/ iv^^p < / IV 



r^llnAI 



It remains to check that the energy contribution related to (p^^^ is of lower order. Indeed, from the 
definition of and ip'^'^'^ it follows that (p*^^^ satisfies 



(1) _ 



in Ztrns 



and 



(0) 



on l9_LZtrns- 
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In particular, at the boundary of the ceh we have 



(1)1 < 



E 



Pi\diPi 



ha 



It then follows by standard elliptic estimates that 



,3.51|l „4 

(1)|2 f r 



on d^^Z . 



r-^llnAI 
h ■ 



This completes the proof of Lemma 3.9 



(3.51) 



□ 



Lemma 3.10 (Energy of a closure domain cell for u^^). Consider a closure cell Z^s with height 
h and maximum needle radius r. Suppose that the needle is slender in the sense ofr<tih. Then 



u 



AF ~AF 



^ ] < erh 



InAI 



+ \r^h. 



(3.52) 



Proof. The estimate for the energy in the closure domain cells proceeds similarly to that for the 

) is give: 

Pcis(a;i) 



transition cells. We decompose ipi — ipf'^ + </?^^', where the function is given by 



ip^°\xi,x±) = 2pcis{xi) diPcis{xi)H(\xj_\ - pc\s{xi))\n- 



\^±\ 



As in the proof of Lemma 3.9 it can be shown that the contribution of the energy related to (ys^^-* 
can be neglected. Hence, in the following, we only give the estimates for the functions ipf'' and (p2- 
A straightforward calculation then yields that we have the following bound on the surface energy 



|VmJ'^^| < erh, 



where we used the assumption r <ti h. Similarly to the arguments in the previous lemma, one can 
also show that 



(0)|2 



< 



I InAI 



and 



(0)|4 < r 



Zeis 



6 
h3- 



Furthermore, from (3.38) and (3.48) we get |9i932| ^ A, and so so 



In view of (3.53) and (3.54), the stray field energy is estimated by 

r J3.33| 



This concludes the proof of Lemma |3.10| 



(3.53) 
(3.54) 

(3.55) 

□ 
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We next give the estimates for u'^^: 



Lemma 3.11 (Energy of a transition cell for u ). Consider a transition cell Ztms with height h 
and maximum needle radius r. Suppose that the needle is slender, in the sense ofr<^h. Then 



Qr-^|lnA| 
h 



(3.56) 



(0). 



-if 



(1) 



Proof. We use the notation of the proof of Lemma 3.9, in particular, we decompose f ~ (p 
with ((3*^"^ being the dominant term. We first note that by construction |V_L(p| < r/h <C 1 (see 
(3.49l), so u"^^ is well-defined in Ztms- Therefore, the surface energy can be estimated as before: 



S,SF< 



IVm'J'^^I < erh. 



Ztr„= 



The estimate for the anisotropy energy of u'^^ follows from 



,SF\2 



Ztr 



llnAI 



h 



(3.57) 



(3.58) 



In order to estimate the stray field energy of u'^^, we note that in view of (3.39 ), we have |v"^^(a;)| < 
|V_L<^(a;)p. Hence 



h^' 



The above estimates together yield (3.56). 



□ 



Lemma 3.12 (Energy of a closure domain cell for u ). Consider a closure cell Z^s with height 
h and maximum needle radius r. Suppose that the needle is slender in the sense ofr<tih. Then 



^|z,Ju^^,v«^] < erh- 



Qr4|lnA| 



+ 7^ + Xr^h. 



Proof. As in Lemma 



3.11 



the function u'^-'^ is well-defined, since |Vj_</ji| < r/h <^ 1. Using the 
slcnderness condition and the fact that by construction |9i(p2| ^ A ^ 1, we again get (3.57). 
Similarly, (3.58) also holds for the anisotropy energy. Finally, in view of the definitions (3.36), 
(3.37), (3.38) and (3.41) and in view of (3.53) and (3.54), the stray field energy is estimated by 



ySF\2 < 



|9l¥'2| 



/l3 



+ Xr'h. 



The above estimates together conclude the proof of Lemma |3.12[ 
We are ready to give the proof of Lemma |3.8[ 



□ 



Proof of Lemma 3.8. In view of (3.46), it is enough to give an estimate in terms of the energy E 
instead of E. As before, consider either a transition or a closure cell Z with dimensions /i, a and the 



maximum needle radius r. The energy of the configurations u , v or u'^ , v'' are estimated 
in the previous Lemmas. In the following, we will write u, v for a configuration representing either 
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^AF^^AF ^SF^^SF^ j^g-j- nionicnt assume that all needles are slender, in the sense of 

r ^ h. In view of Lemmas 3.9 3.12} we then have for any transition cell Ztms and any closure 
domain cell Zeis with the above dimensions the following estimate: 



7r'*|lnA| 

h ^ ¥ 



and i?|z^jju,v] < er/i- 



jr^ I In A I r 



+ -r^ + Xr'h. (3.59) 



Balancing the first two terms in the right-hand sides of (3.591 yields the optimal height for both 
transition and closure domain cells as a function of the maximum needle radius r: 



72 £ 2 7-b|l];iA|2. 



(3.60) 



With this choice of h for each cell, the estimates in (3.59) turn into 

„6 _ _ 



|lnA|5 



r 



and E\z^i^[\i,w] < j^e^r^\liiX\^ + + Xr^h. (3.61) 



Given any initial height hi for the first generation of cells, we use (3.60) to correspondingly choose 
the width of the first generation of cells. The width of the following generation of cells is inductively 
defined by (3.22) and (3.60). We terminate the algorithm after M generations of cells as soon as 
closure domain cells are not too expensive in the sense of 

74/1 In A I 



XrlihM 



< 



(3.62) 



We choose the initial height hi = hi{M) such that (3.23 ) is satisfied, i.e. such that the cells exactly 
cover Q. Since in view of (3.22) and (3.60), hj is a geometric sum, as expected we must have /ii ~ 1 
independently of M. In view of (3.60) and (3.62), we then get the following estimate for the needle 
radius for the first and last generation of cells 



ri 



£3 



73 1 In A| 3 



and 



rM 



Note that the termination criterion (3.62) is equivalent to rfj/hj^j 
(3.22) and (3.60) Vj/hj is monotonically increasing in j, we get 



< 



A 



7|lnA|- 



< 



A 



7|lnA| 

Let us assume for the moment that 



< 



7r4|lnA| 



for all < j < M. 



for all < j < M. 



In this case, in view of (3.61 1, the total energy is estimated by 
1 ~, ^ 1 ^ 1 1 ^ ^ ^ 



InAh 



7363 A| lnA|: 



(3.63) 

Since in view of 
(3.64) 

(3.65) 



25 



In order to complete the proof, it remains to check the following three consistency criteria. 



We first need to verify (3.65). Indeed, (3.651 follows from (3.64) and our second assumption on A 



Secondly, we need to check that our algorithm allows for at least one generation of cells, i.e. Af > 1. 
In order to see this, we note that M > 1 is equivalent to 3*^ > 1. By our first assumption on A, we 
have 



A 



73 £3 1 In A| : 



> 



A 



7 3 e 3 1 In £ I 



> 1. 



Finally, we need to check that the cells are indeed slender in the sense of rj <^ hj. Indeed, this 
follows from the second of (3.64), together with the second assumption on A. This concludes the 
proof of Lemma |3.8| □ 



3.5 Constructions for the full energy 



In this section, we give the proof of Theorem 3.2 In order to do so, it remains to give a diffuse 



interface version of the upper constructions in Section [3^3] 

We first consider the case of a hard material, i.e. Q ^ 1. Hence, we will construct a diffuse interface 
version u'^^ of the magnetization u^^. It is enough to show the construction for a single needle. For 
simplicity, consider a closure cell Z^is with height h and width a. We recall that for sharp interfaces, 



we defined u'^^ by (3.32), where the shape of the needle is described by the characteristic function 



X in (3.30). Now, we define the transition layer 

S = {{xi,x±_) e Zeis ■■ \ \x±\- pcis{xi) - d{xi)\ < w{xi)}, (3.66) 

where the functions w{xi) and d{xi) are the thickness and the displacement of the diffuse interface, 
respectively, given by 

w{xi) := ^ Pcis{xi), d{xi) := -pcis(a;i) (l - ^1- 3^1^^- ) • (3.67) 

We recall that in the present units the quantity eQ^^ is just the typical width of Bloch walls, as 
described in Section [2] We also note that by our assumptions eQ~^ <C r and, therefore, |(i(a;i)| <C 
w{xi) <^ PcIb{xi), for all xi € (0,/i) and all cells. Indeed, by our assumptions 7 > X^\lnX\~^ » A. 
Therefore, by ( |3.63[ ) the inequality holds for the closure cell, in which r is the smallest. 

Let the "mollification of the Heaviside function" H € W^^°°{M.) be given by H{t) — in (—00, — 1], 
H{ t) = ^ {t+ 1) e [-1,1] and H{t) = 1 in [l,oo). Furthermore let HR{t) := H{t/R). Analogously 



to (3.30), we define 

X{xi,xj_) H^^^^){pcisixi) + d{xi) - \x±\) for x G Zeis, 

We then define u"^^ as follows. First we set 

iif^ = A — 2x{x) in Z^s and uj^ := outside of S. 

Inside the transition layer, we define u^^, such that |u^''^| ensures that u £ A, and that the vectors 
ei, x±, and form a right-handed triplet. A straightforward calculation shows that our choice 
of d{xi) ensures charge neutrality on every slice, i.e. Jj^ uf^^ {xi, x ^) dx ~ for all xi G [0,h]. 
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In the above construction the exchange and anisotropy energy are supported only in the transition 
layer S. One easily gets that 



4Q 



Q\u_ 



< 



e'^Pcisjxi) 



Qpcisixi)w{xi) dxi 



erh. 



In estimating the stray field v of u , we can follow the same arguments as for u , modifying 
the definition of ip^^^ to be the radially-symmetric potential due to charges V • u^''^ rather than 
V • u'^-^. It is easy to see that all the estimates remain unchanged. Comparing the above estimates 
with those in (3.52), it follows that the localized diffuse interface energy £[u'^^, v'^-'^] for a single 
cell, based on the construction u'^^ is not larger in terms of scaling than the loca lized sharp interface 
energy E[\j/^^ of the optimal sharp interface construction in Section 3.3 



The function u'^^ can be defined throughout by ap plying the above construction to every needle 
in the self-similar geometry described in Section [3^ The corresponding estimate for u"^^ follows. 
Lastly, we note that the construction in the case of a soft material {Q < 1) proceeds analogously. 



This concludes the proof of Theorem 3.2 



4 Reduced energy 

The analysis of Section [3] provides the scaling of the minimal energy in the limit of thick samples 
or, correspondingly, when e — > 0. While the analysis does not require any assumptions about the 
minimizers (the analysis performed is ansatz-free), the results obtained give only a rough idea about 
the structure of the minimizers. It seems natural to expect that the minimizers should look like 
the trial functions used in the construction of the upper bounds. Yet, the precise shape of the 
domains, as well as the precise constants in the asymptotic behavior of the minimal energy cannot 
be captured by the analysis above, since it does not address the leading order constant in the scaling 
of the energy. 

In this section, under the assumption that the magnetization is mostly aligned with the easy axis 
and that the geometry of the minimizers is slender, we derive a reduced sharp interface energy 
which should provide the leading order behavior of energy for e <^ 1. Our aim is to reduce energy 
minimization of i? to a two-step process. In the first step, we fix the "shape" of the magnetic 
domains and construct the energy-minimizing configuration of the magnetization away from the 
domain walls. In the second step, we minimize the obtained, reduced energy, which depends only 
on that shape. We note that the heuristic idea of computing the combined contribution of the 
magnetostatic and anisotropy energies away from the domain walls has been known as the /x*- 



method in the physics literature 16 33 . Below, we assign a precise mathematical meaning to this 
idea and provide its rigorous justification under specific assumptions. Finally, let us also point out 
that while in this paper we are interested in the case of applied field near saturation (A ^ 1), we 
expect that the obtained reduced energy to be valid independently of the applied field, even in zero 
applied field (A = 1). 

We first introduce the characteristic function x representing the shape of the domains where the 
magnetization vector is not aligned with the external field, 

X{x) = 1 if mi(a::)<0, x{x) = if mi{x)>Q. (4.1) 
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The reduced energy which is derived in this section is then given by 

Eo[x] = 2e f |Vx| +7 / di{Xxn - 2x)(-AQ^)9i(Axa - 2x), (4.2) 

where the operator Ag is defined by 

Aq = (xo + 7(1 - Xo)) A_L + jdf, 



with constant 7 defined earher in (3.1 ). The admissible class of functions x for ^0 is 
Ao = {xe BV(R X T; {0, 1}) : x = in (M x T)\n}. 

Note that Aq « Aj^ in Q, when acting on functions that vary slowly in the easy direction, compared 
to the directions normal to the easy axis. There is a second equivalent formulation for (4.2 ). Indeed, 
a straightforward calculation yields that (4.2) can be written as 

Eo[x] = X^i^ + 2e f |Vx| - 4A / x + 47 / d,x{-^Q')diX- (4.3) 

Notice that by lower semicontinuity and coercivity, the minimum of the energy of Eq is attained in 
Aq (see also [6j Theorem 1.2]). 

We first note that up to the leading order constant, the reduced energy Eq has the same scaling of 
minimal energy as E, i.e., for A < 7^| In Ap and £ sufficiently large we have 



4 inf Eq[x] - min|A^73^^lA|lnA|3| 



It can be checked that this result follows by a slight modification of the proof of Theorem 3.1 
(replacing ui with Axo — 2x). 

In order to show the asymptotic equivalence of the minimum energies for E and Eq, including the 
leading order constant, we need to make an assumption on the magnetization m. If we assume that 
the magnetization vector m does not deviate strongly from the easy axis throughout the sample 
and, furthermore, that the geometry of the magnetization configuration is slender, we can show 
that the minimal energies for Eq and E essentially agree to the leading order: 



Theorem 4.1. Let £ <C 1 and j ^ S, where < (5 <C 1 is the same as in (2.16). Then 
(i) For every \i ^ A satisfying |uj^| < 6 there exists x G -Ao, such that 

E[u\ > {1 - S-2)Eo[x]- 
(a) For every x G -^o, for which the solution (p of 

AQ(^ = 79i(Axn-2x). (4.4) 
satisfies |Vi^| < Q6, there exists u Cz A, such that 

Eo[x] > (l-(55)^[u]. 
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Theorem (4.1) is proved at the end of this section via Propositions 4.2 and 4.3 Before giving 
the proofs for Propositions 4.2 and 4.3 it is instructive to present a formal derivation of Eq (see 
also 33 ). 

We note that both anisotropy and the external field favor alignment of the magnetization with the 
easy axis. We hence expect that m ^ ±ei in and should, therefore, have wi ~ A — 2% there. In 
view of |m| = xn, this motivates to write u in in the form 



A - 1 + (1 - 2x)v/l - |u±P ei + = (A-2x)ei+Ui + 0(|ui| 



For |mj^| = |uj^| <C 1, to the leading order (2.161 then turns formally into 

E[u]c,2e f |Vx|+Q / |u^P+ / |vp. 



Here, the function v satisfies 
V = 



A(p = ai(Axn-2x)+V_L-u_L 



X T. 



(4.5) 



(4.6) 



Following our approach, we minimize (4.5) in two steps. First we take the minimum with respect 
to u_L with X fixed. In the second step we minimize the result with respect to all admissible 
characteristic functions x- It is not difficult to see (see the proof of Proposition 4.2 for details) that 
for fixed x the sum of the last two terms in (4.5 ) is minimized when 



Substituting this relation into (4.6), we find that 

dl^ + A^<f^ di{Xxn-2x) + y±-n^_ = di{Xxn - ^x) ^ Q^'xnA^^, 



which is precisely (4.4). Finally, substituting the expression for in (4.4) into (4.5), it then follows 
that £'[u] ~ -Eoixli where the reduced energy Eq is given by (4.2). Here we took into account that 
Aq is an invertible operator. Thus, to the leading order the energy of minimizers of Eq should 
coincide with that of E. 

Observe that the reduced energy £^o[x] j^st derived has a form very similar to that of the sharp 
interface energy E in the case of infinite anisotropy, Q — oo. Indeed, in the latter case the mag- 
netization vector m is restricted to take only two values in f2: m = ±ei. For such magnetization 
configurations the sharp interface energy (2.16) turns into 



Jn JRxT 



which coincides with (4.2 1 for 7 = 1, since in this case ui — Xx^ 
that for slender magnetization configurations we have — A^^ ~ 



2x. On the other hand, noting 
-A]]^, one should expect 



E 



AF 



2e 



\^±X\ + 



9i(Axo-2x)(-Ali)9i(Axo-2x), 



Eo^2e [ |Vixl+7 / M^Xi. 
Jn Jrxt 



2x)(-Ali)9i(Axo-2x), 



(4.7) 
(4.8) 
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where x in (4.7 1 and (4. 



is given by (4.1). It is easy to see that up to a muhipUcative factor 

Similarly, 



the expression in (4.8) coincides with that in (4.7 1 after rescaling X2 and 2:3 with 7 



the energy per unit area, according to (4.8), is 73 times the expression in (4.7). Thus, the two 



energies approximately agree with each other for Q 3> 1, and the minimal energy per unit area is 



smaller by a factor Qa for Q ^ 1 (see (3.1)). The latter is due to the fact that in this case the 
effective magnetostatic interaction is weakened by a factor of Q, since the stray field is shielded by 
small deviations of m from the easy direction, creating magnetic counter-charges. These arguments 
provide a physical explanation of the apparently surprising fact that the scalings of the energy of 
non-trivial minimizers both in the case of hard and soft materials agree up to a factor of (see 
also 1^^). 

We now give a rigorous derivation of the relationship between the reduced energy Eq and the sharp 
interface energy E under a few assumptions which appear quite natural physically. Theorem |4.1| 
is an immediate consequence of the two propositions that follow (for the lower and upper bounds, 
respectively). We begin with the analysis of the lower bound for E[xi\ in terms of £^o[x]- 



Proposition 4.2. Let 7 ^ (5, where < (5 ^ 1 zs the same as in (2.16), let u d A, and let 
|u_l| < S. Then E[u] > (1 ~ S^)Eo[x\, where x is given by 

u(0) 



Proof. Let us write u 



u 



(0) = 



+ u^^^ and (p — i^'"'' -I- ip''^\ where ip is defined in (2.17), and 



= X^^{(A-2x)el-Q-lVi^W} 



(4.9) 



This is always possible, sinc e (p^^^^ is uniquely solvable in terms of x- Indeed, eliminating u^^-' in 
the second equation in (4.9) via the first equation in (4.9), one immediately sees that (^^°^ solves 



the same equation as (f in (4.4). Then, after a straightforward computation, using (2.17) and (4.9) 
the bulk part of energy E'buikiu] = Q Jq |u_lP -I- /jj^t I^*^!^ written as 



E^ 



bulk 



U] =Q 



u 



(1) 



(0) 



(1) 



Q 



.(1) 



(0) 



2^1 v^^"^ u^'A + 



Using Cauchy-Schwarz inequality, we obtain 



-Bbuik[u] > / (l + 0"'xf2)|Vj.(p 



(0)|2 



\diip 



(0)|2 



\diip 



(0)|2 



(0) 



(1) 



,(1)|2 



(4.10) 



ei, where rn^ is the first component of the vector m'^°) — u^^^^ -I- (1 — A)ei, does not exceed 



Now, recalling (2.15 ) and (4.1 ), one can see that in O the angle between the vector m and the vector 



^. Note that by (4.9) we have = 1. Therefore, since |m| = 1 as well, the magnitude of the 

projection of the vector u*^^-' = m— m^"' onto ei does not exceed 1, and u^^'ei points in the direction 



opposite to mf'^'ei. From this we conclude that \u\''^\ = 1 — ^/l — \u±\^ < \u±\'^. Therefore, by 
assumption we have {u^i'^l < 6\u±\. On the other hand, as can be easily seen, the non-local part of 
Eq equals the sum of the first two terms in (4.10). Furthermore, by Cauchy-Schwarz inequality 



diip 



(0) 



,(1) 



< 2 Eo[x]S^ / |ui 



< 2SQ~-2Eo[x]-^Eo[u]-2. 
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Combining the above estimates with the fact that 

/ \Vut\ = 2(1-6') [ |Vx|, 
Jn Jn 

we arrive at 

E[u]>2eil-d') f |Vx|+£;buik[u] >{1-S')Eo[x]~2SQ-^E-2[u]e'^[x], 
Jn 

which yields the statement. 

We next give the proof for the upper bound for E[u] in terms of -EqM- 



□ 



Proposition 4.3. Let x G -^o, let |Vj^i^| < Q5, where (f is given by (4.4), for some < (5 ^ 7. 
Then there exists u G ^ with \uj_\ < 6, such that Eq[x\ > {1 — S^)E[u]. 

Proof. Set ^ ^ and u = u(o) + u(i), where u'^) = xni{>^ - 2x)ei - Q^^V^l'^^''^), and u^^) = 
u^^'^-'ei ensures that ju^''^ + u'^-'^^ + (1 — A)xnei| = xn- Then (p*^") solves Aip^^') — V • u^^^^, and by 
assumption lu^^"*] < Q^'^ l Vj_yj^°-' p. The proof is then obtained by retracing the calculation in the 
proof of Proposition 4.2 noting that in this case u^' = 0. We will only need one extra estimate 
for ip'^^^ solving Aip'^ ' — V ■ u'--^\ 

Integrating by parts and applying Cauchy-Schwarz inequality, we obtain 



<T Jmxt Jrxt Jn \Jn Jn 



Squaring both sides and using the above estimate for 



,(1) 



(1) 



(0) 



, we find that 



< S'Q-^Eoix], 



which completes the proof. 



□ 



5 Transition to non-trivial minimizers 



As we showed in Theorem |3.2[ there is a change in the scaling behavior of the minimum energy 
due to appearance of non-trivial minimizers of £ at A '-^ 73 ea | lne| 3 . In this section we analyze 
the nature of this transition in more detail. Specifically, we are interested in locating the precise 
critical value of A (corresponding to the critical applied field away from saturation) at which this 
transition occurs. We also address the structure of the domain patterns near the transition point. 

As a first step, based on an asymptotic study of the reduced energy Eq, we derive an even further 
reduced energy £'00 ■ We will formally show that near the transition point, in rescaled variables, it 
is appropriate to consider the energy 



Eoo[A] 




A- 



(5.1) 
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where the set of admissible functions is 

A = {Ae i?^((0,/i)) : A > 0, A(0) = = 0} 

The function A{^) is simply the rescaled area of the cross-section of a single needle as a function of 
the rescaled coordinate along the needle. The single parameter h can be understood as a measure 
for the effective thickness of the plate. It is defined by 

At 

(5.2) 



£72 In 2 (732 



A detailed derivation of the reduced energy Eqq and the precise definition of the quantities A and 
^ in terms of the original quantities is given in the next subsection. 

The advantage of the energy i^oo is that it can be explicitly minimized and its minimizers can be 
explicitly computed. We identify two critical values, denoted as /ij and h\, with < h\ for the 
thickness of the rescaled slab, at which transitions in the qualitative behavior of the critical points 
of -Eqo occur. Basically, the result is that the uniform state is the unique global minimizer whenever 
the effective thickness h of the sample satisfies h <h\. Furthermore, the uniform state is even the 
unique critical point as long as h < h^. The precise statement is the following: 

Theorem 5.1. Let =^ TTy/2 and let > Iiq be the unique solution of the system of equations 

F{pm)=0 and G{pm) = hl (5.3) 



where the functions G and F are defined in (5.8) and (5.9), respectively. Depending on the value 
of h, we then have 

1. If h < /iq, then A = is the unique global minimizer of Eqq, and there are no other critical 
points of Eqo . 

2. If h^ < h < hi, then A — is the unique global minimizer of Eqq, but there exist non-trivial 
critical points of Eqq . 

3. If h = hi, then there are two global minimizers, given by A — and by the unique positive 
solution of (5.5) vanishing at the endpoints. 



4. If h > hi, then the unique positive solution of [5.5) vanishing at the endpoints is the unique 
minimizer of Eqq . 

Furthermore, the unique positive critical point Ah of Eqq obeys Ah{S,) ~ C"*/"^ for h — and 
Ah{0-Uorh>h*Q. 

Thus, the transition from trivial to non-trivial minimizers of Eqq occurs precisely when h — hi. 
Numerically, the critical values of the parameter h are 

h*Q w 4.443, hi w 6.113. 



In terms of the original energy £, the statement of Theorem |5.1| has the following interpretation. 
In view of (5.2), the critical values /ij and hi of h define the respective critical values of A: 
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with the meaning that if one chooses A = Aes | In e| 5 ^ then the global minimizer of £ will be trivial 
when A < A^, and non-trivial when A > A^, for fixed A and 7 as e — > 0. Similarly, the trivial 
minimizer is expected to be the unique critical point of £ for A < Aq when e — >■ 0. Thus, the 
transition to non-trivial minimizers for e ^ 1 is expected to occur at A = A^; — A^^e^llneja. We 
would similarly expect the transition to non-trivial minimizers to occur at this value of \*i in all 
the energies: i?o, and S. We note that A^ can be rigorously shown to give the asymptotic upper 
bound for the critical value of A at which non-trivial minimizers emerge by constructing suitable 
trial functions out of the non-trivial minimizers of Eqq. However, since the arguments in this section 
are based on certain assumptions on the geometry of the magnetic domains, the arguments are not 
rigorous in terms of a lower bound for A^^ . 



5.1 Isolated needles 



In this section, we present a formal asymptotic derivation of E'oo in (5.1). At the onset of the 



transition from uniform magnetization to a patterned state as the applied field is reduced it seems 
natural to expect the appearance of thin slender needle-shaped domains of magnetization opposing 
the applied field. Under this assumption, it is possible to further reduce the energy Eq to obtain 
the precise information about the shape of these domains. 



The starting point of the analysis in this section is the reduced energy Eq in the form of (4.3). We 



are interested in the magnetization configuration in the form of a single needle. More precisely, 
we assume that the configuration consists of a a single needle-shaped domain in a sufficiently large 
sample (i.e. ^ ^ 1). In particular, suppx looks like the characteristic function of a prolate ellipsoid 
of radius tq <C 1, extending across 17 in the direction of the easy axis. The crucial observation for 
the analysis in this section is that due to the slender geometry, the dominant stray field interaction 
is restricted to slices normal to the easy axis. This interaction, however, is logarithmic and hence 
does not see the precise shape of the magnetic domains (a similar phenomenon occurs in a related 
model 26 ). In fact, for a domain pattern described above, to the leading order in ro ^ 1 the 
Green's function Gq of the operator — Ag can be approximated by Gg given by 



Go(r) 



27r 



Indeed, Go gives the leading order behavior of the Green's function for the operator — Aq 
Then the non-local term in the definition of E^ may be written as 

JRxT Jo Jt Jt 



^A, 



|hiro 
2 



This motivates to define the function ^ : [0, 1] — >■ M by A{xi) — /rpx(2^ij')j denoting the cross- 
sectional area of the needle in the slice at xi. Then (4.3) turns into 



i?o[x] « A^^^ + 4AA + ^l\}p\\Af^ dxu 



33 



where we again have used slenderness of the needle in the sense of |Vx| ~ |Vj^x|. Minimizing 
the interfacial contribution of the energy at each xi for fixed cross-sectional area then leads to the 
following expected behaviour for minimizers: 



inf Eq [x] 

X 



inf£:oo[^], 



where -Eqo is given by 



Em[A] = / \Ae^/^TA'^ - AX A 



27|lnro| 



Af ]dxi. 



(5.4) 



Note that according to ( 3.63 ) , at the transition, where A 
and, hence, |lnro| ~ |ln(7/e)3 



73£3 1 lne| 3 , we expect ro ^ £37 3|ln£| : 
Using these two scalings and by introducing the rescaled variables 



.- _ A2 
A--^A, 



^ = hxi, 



^oa[A] = ^Eoo[A], 



we obtain the energy Eqo in (5.1 ). Finally, since ^(0) > or A{h) > imply that there is a charge 
layer at the plate's surface, causing a lot of stray field energy, to the leading order the minimizers 
of Eqo are expected to satisfy A{0) = A{h) = 0. 



5.2 Needle shapes and critical fields 



In this section, we investigate the minimizers of £^00 and give the proof of Theorem 5.1 
assume there exists a positive critical point Aii{£,) S HQ{{0,h)) of Eqq 



I.e. 



We first 
Ah>0 in {0,h). By 



standard ODE theory, Ah £ C^{{0, h)) satisfies the Euler-Lagrange equation for (5.1), i.e. 
1 d^Ah 



7T de 



= -1 



2VAh 



Note that (5.5) admits a first integral 

2 



1 fdA, 



2tt V 



Ah - \f^h = ttC 



for ^ e (0,/i). 



for i e (0,/i), 



(5.5) 



(5.6) 



where C G R is an arbitrary constant. Evaluating (5.6) at = and in view of AhifS) 
C > 0. Evaluating (5.6) at the maximum point of A/j, yields maxA/j 
that the solution of 



> TT. 



0, we have 
We furthermore note 



5.6) is monotone in (0, ^) and takes its maximum at ^ 



■| . It is convenient to 

introduce the rescaled needle radius p(^) = (^ft(^)/7r)5 where in view of the above the maximum 
Pm = p{h/2) satisfies pm > 1. Integrating (5.6) over (0, |), a straightforward calculation yields 

h 
2 



V2 



pdp 



^C + p-p^ 



2pni) + 2^J Pra{Pra - 1)) • 



(5.7) 



Here we used the fact that C — Pm{Pr. 

particular, for any given /i > a positive critical point of £^oo exists if and only if h 
some Pm > 1, where 



1), which follows by evaluating (5.6)at^ = ft,/2. In 

G{pm) for 



G(p,„) := ^{^QC ^(1 - 2p,„) + 2\J Pra{Pm - 1)) 



(5.. 
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Figure 6: a) Parametric dependence of needle height on the radius, obtained from (5.8). b) Needle 



energy as a function of radius, obtained from (5.9) 



Differentiating (5.8), one gets that '^'^Jf 
increasing for p. 



^ _ 4^2p„(p„-l) 



2p > 0, i.e, G{pm) is strictly monotonically 

> 1, with G(l) = ■K\/2 and G{pm) ^ oo as pm oo. In particular, a non-trivial 



critical point of i?oo exists, if and only \i h > Hq := n\/2. Since for every local minimizer the 
Euler-Lagrange equation holds in every interval of positivity, this also shows that for every h < h^, 
the only critical point of i?oo in ^o((0; ^)) is ^ = 0. This completes the proof about the existence 



or non-existence of non-trivial critical points in (i) and (ii) in Theorem 5.1 



Let us now consider the global minimizers of i^ooi which exist, in view of coercivity and lower- 
semicontinuity of i?oo for all /i > 0. We first calculate the energy of the needle profile calculated in 
the first part of the proof. For this, we define 



An explicit computation then yields 



dAh 



- irhC. 



F{pm) = -^(3 - 4p„(p„ - l))v/pm(Pm - 1) + ^ SCC ^(l-2p,„). 

V 18 V 8 



(5.9) 



rff (p„) _ 47r^2p-^(p„-l)3 



< 



Once again, differentiating this function with respect to pm, one gets — ^ — 2p -i 

0, so that F{p,n) is strictly monotonically decreasing for p„ > 1, with F{1) = 7r^/(2-\/2) and 
F — > — oo as ^ — > oo. By monotonicity of G it then follows that i5oo[^ft] is strictly decreasing in h 



for h>h*Q. In particular, for h\ > defined by ( |5.3[ ), we have ii^oo[^/i] > for h e {hQ.hl) and 
£'oo[^ft,] < for all h € {hi, oo) In particular, for h G (0, /ij), the only global minimizer of £"00 is 
given by ^ = 0. The dependences of h and i?oo[^/i] on p„i are shown in Fig. |6] 

We now claim that for h > h\ the global minimizer A is unique and is given hy A = Ah. Indeed, 
we first note that A is not equal to 0. In view of the above estimates on _Eoo[^?i]) every interval of 
positivity contains a point x* with A{x*) — tt. By strict monotonicity of i?oo[^/i] as a function of 
ft., it is furthermore clear that ^ = and £, = h are boundary points of the intervals of positivity of 
A. Suppose that A = Q on I = [^1,^2] CC (0, /i) where ^1 < ^2- By the above reasoning it follows 
that there exist points G (0,6) and £,2 G (67^) such that A{£,{) = ^(^2) = ^- K follows that A 
defined by A := tt in (^^,^2) and A:— A outside of (^11^2) -^las lower energy than A contradicting 
the assumption that A is a minimizer. This shows that A > in (0, /i) and hence A = A-^. This 



completes the proof of (i)-(iv) in Theorem 5.1 
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Figure 7: The needle shape at h = hi, obtained from (5.5) 



In view of (5.6 1 and (5.7), an explicit calculation yields the following parametric equation for the 
needle profile for < ^ < |: 

^ = (2v/(p™ - l)p,n + tan-1 ((2p - l){2^ip,n ~ p){pm +P- 1))"') 

see Fig. [Tj In particular, one easily checks that for < ^ <C 1 the behavior of the radius p of the 
needle near the tip is given by 

p(^) ^ ^5 for/i>/i* and p{C) ~ for /i = /i*. (5.10) 

This concludes the proof of Theorem |5.1[ 
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